Why spin ice obeys the ice rules 



in 
o 
o 

(N 

Oh- 

<D 

in 



o 



CO 



I 

C 

O 

o 



> 

m 

O 
in 
o 



-a 
c 

o 
o 



X 
S3 



S. V. Isakov, 1 R. Moessner, 2 and S. L. Sondhi 3 

1 Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 
2 Laboratoire de Physique Theorique de I'Ecole Normale Superieure, CNRS-UMR8549, Paris, France and 
3 Department of Physics, Princeton University, Princeton, NJ, USA 

The low temperature entropy of the the spin ice compounds, such as Ho2Ti207 and Dy2Ti207, is 
well described by the nearest-neighbor antiferromagnetic Ising model on the pyrochlore lattice, i.e. 
by the "ice rules". This is surprising since the dominant coupling between the spins is their long 
ranged dipole interaction. We show that this phenomenon can be understood rather elegantly: one 
can construct a model dipole interaction, by adding terms of shorter range, which yields precisely 
the same ground states, and hence T — entropy, as the nearest neighbor interaction. A treatment 
of the small difference between the model and true dipole interactions reproduces the numerical 
work by Gingras et al. in detail. We are also led to a more general concept of projective equivalence 
between interactions. 

PACS numbers: 75.10.Hk, 75.50.Ee, 75.40.Cx 



Introduction: In 1956, Anderson observed that an 
Ising antiferromagnet on the pyrochlore lattice would 
exhibit a macroscopic ground state entropy equivalent 
to that of water ice 0]. Almost four decades later, in 
1997, such "spin ice" behavior was experimentally dis- 
covered. Harris et al. noticed that the pyrochlore com- 
pound Ho2Ti207, where only the Holmium ions are mag- 
netic, failed to display any signs of ordering down to be- 
low the ferromagnetic Curie temperature deduced from 
its high temperature susceptibility. They proposed that 
this was due to the presence of an easy-axis anisotropy, 
which leads to the ferromagnet effectively becoming an 
Ising pseudospin antiferromagnet at temperatures well 
below the anisotropy strength This interpretation 
was strongly supported by a remarkable experiment, in 
which Ramirez et al. showed that the related compound 
Dy2Ti2C>7 displayed a residual entropy, So, at low tem- 
peratures, which was in excellent agreement with the 
Pauling estimate for water ice, So « (1/2) ln(3/2) 0. 

However, it was quickly pointed out that the dominant 
interaction in these compounds is dipolar due to the large 
moments on Ho 3+ and Dy 3+ , /i w 10/is, where hb is the 
Bohr magneton Yet, the dipole interaction is long 
ranged, decaying with separation, r, as 1/r 3 . Why, then, 
is the entropy of the nearest-neighbor model stable to the 
inclusion of the rest of this interaction 5, 61? 



That this is the case has been checked in detail by ex- 
plicit simulation. Further, two qualitative observations 
have been made in the literature. First, the anisotropic 
nature of the dipole interaction precludes any obvious 
ordering instability stemming from its long range. It is 
probably fair to say that it is now agreed that instead, 
it leads to a weak ordering instability away from q = 
HUSH], albeit one which is not observed || . Sec- 
ond, explicit evaluation of the Fourier transform of the 
easy-axis projected dipole interaction to large (and, us- 
ing Ewald summation, infinite) distances yielded a sur- 
prise: the infinite distance result resembled the nearest- 



neighbor interaction more than an interaction truncated 
at, say, 10 or so nearest neighbors; this 'self-screening' 
(Gingras) indicates that there is something special about 
the dipole interaction 0, 0] . 

In this paper we show what it is that is special about 
spin ice and dipole interactions. The main insight comes 
from recent progress in understanding the correlations 
dictated by the ice rulejll El El Q: it was shown 
that the pyrochlore Ising model dynamically acquires a 
gauge structure at T — 0, which manifests itself in the 
emergence of dipolar correlations as the ice rules are en- 
forced. The crucial observation is that this gauge struc- 
ture is (formally, not in origin) exactly the same as that of 
'ordinary' magnetostatics; this latter, of course, is what 
determines the form of the dipole interactions. 

Most crisply, we use this observation to show that 
there exists a slightly modified "model" dipole interac- 
tion, which differs from the physical interaction by terms 
that fall off faster and are small in magnitude, whose 
ground states are identically the same as those of the 
nearest neighbor interaction. This accounts for their 
identical low temperature entropy. 

The mathematical implementation of this insight re- 
volves around relating two quantities - dipole interac- 
tions, and the operator enforcing the ice rules - to the 
same projector, V . Finally, this leads us to a generaliz- 
able equivalence between a long ranged interaction and 
one of shorter range projected onto its low energy mani- 
fold; this we term "projective equivalence" . 

In the balance of the paper we flesh out these state- 
ments. We begin by recalling the energetics of spin ice. 
The ground state degeneracy of the nearest-neighbor an- 
tiferromagnet manifests itself in a pair of flat bands. We 
show next that a model dipole interaction possessing the 
correct long distance form, leads to exactly the same 
pair of flat bands and thence to the same ground state 
manifold, even with arbitrary admixture of the nearest- 
neighbor interaction (as long as not so overwhelmingly 
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antiferromagnetic to invert the ordering of the bands). 
The deviation of the model interaction from the true 
dipole interaction is seen to vanish as an integrable r~ 5 . 
Treating this by elementary means we are able to account 
for the weak residual dispersion and ordering tendency 
found in previous work. As a byproduct, this analysis 
makes clear that at higher temperatures, the correlations 
arising from the dipole interaction are not the same as 
those stemming from the ice rules alone. We close with 
some comments on projective equivalence. 
Hamiltonian: The dominant term in the spin ice Hamil- 
tonian is a strong easy easy axis anisotropy, which al- 
lows us to transform the vector spins S to Ising pseu- 
dospins <7 = ±1 via Si a — ai a e a (no sum over a). 
e a denote the local easy axes of the pyrochlore lattice, 
which consists of corner-sharing tetrahedra. It can be 
thought of as a face-centred cubic lattice with a four 
site basis: ei = (— 1, — 1, — e 2 = (1, — 1, l)/\/3, 

e 3 = (1, 1, —1)/ a/3, e4 = (— 1, 1, l)/\/3; i is a unit cell 
index, and a is the sublattice index. 

The dipolar spin ice model contains two terms. First, a 
nearest-neighbour exchange, of strength J, to which 
we add a constant so that its ground states have zero 
energy; these ground states obey the ice rules or equiva- 
lently the constraint that the total pseudospin of each 
tetrahedron vanish. And second, the dipole interac- 
tion T>, which is summed over all pairs of sites, with 
D = /io/i 2 /(47rrj^ n ), /i is the magnetic moment of the 
spins, r nn is the nearest-neighbor distance, and Tiajp 
is the vector separating spins Si a and Sjp. Defining 
^ = £ia,i0'%a,j/3 (J ia (T i/3> tne Hamiltonian is: 

ti = ^ &ia [jJiajfS + Dr^jjpiajp] (Tjp, (1) 
pairs 

_ e a ■ ep _ 3(e a ■ r iatj p)(ep ■ r iaJI3 ) 

u i<*,3P ~ , ,3 , ,5 W 

Spectrum and eigenvectors of H: In this section, we 
establish the connection of the ice rules with the dipole 
interaction; this we do by discussing the spectrum and 
eigenvectors of the Hamiltonian matrix TL. These, specif- 
ically two dispersionless bands, will form the basis of our 
discussion of the zero temperature entropy and zero and 
finite temperature correlations below. 

First, consider the spectrum of the adjacency matrix 
Jia,j0- It is well know that this has a pair of degener- 
ate flat bands and two dispersive bands (one of them is 
gapless at zero wavevector), see F ig. [H The diagonaliz- 
ing transformation is given in Ref. [l2|. We can write the 
result as the schematic decomposition, 
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.7 = ]r^(q)Mq))Mq)| (3) 
where for /i = 1, 2, e M (q) = 0. 
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FIG. 1: (color online). Mode spectrum of the nearest- 
neighbor interaction matrix, J, in the [hhl] plane fl^l . q is 
in units of 2n. Note that here and on the following plots we 
plot the eigenvalues, Ui, with minus sign. 

As shown below, at T = 0, the physics is determined 
solely by the modes in the zero-energy flat bands; it will 
therefore be unchanged for a family of interaction matri- 
ces with any other choice of e 3i 4 > (but keeping the 
vectors v unchanged). The family member with e 3j 4 = 1 
is then a simple projector V = V 2 = Ylu=3 I^m (l)) (^m (l) I 
(Fig.©. 

This is a l ongi tudinal projector; for formal details, see 
Refs. [H E Il4l Il5l but in a nutshell, what happens is 
this. The ice rules X^etet Gi are equivalent to V • S = 0. 
Here, the spins are thought of link variables ('lattice 
fluxes') on the diamond lattice, which is dual to the py- 
rochlore, and V- is the appropriate lattice divergence. 
This equation encodes the statement that longitudinal 
modes, those with V • S ^ 0, cost energy, whereas trans- 
verse ones do not. The energetically enforced constraint 
V • S = is at the origin of the 'emergent gauge struc- 
ture' of spin ice, as it can be resolved by transforming to 
a gauge field (vector potential), S = V x A. 

What are the matrix elements of V in real space? 
Working backwards from the form of the spin ice cor- 
relations obtained in [l2l H3I ITil Il5| . one can read off 
that, asymptotically at large distances, V oc V . More 
precisly, defining a correction term A through 

87T 

Diajp = —ViaJ/S + &ia,j0 j (4) 

analyticity and symmetry considerations give Ai a jp ~ 
0(r- a %). 

Thus, an interaction V (i) has the same long-distance 
form as T> and (ii) has the same eigenvectors and same 
ground-state manifold as J . With this in hand, we can 
now understand all the important qualitative features of 
the spectrum of Ti. analytically and gain a quantitative 
understanding with computations that do not require us 
explicitly to treat the conditional convergence of dipole 
sums. (Note that we again add an overall constant to the 
energy so that the flat bands occur at zero energy). 
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Starting with the model dipole interaction V (Fig. [21 
top panel), it is trivial to add in the superexchange as J 
and V have the same eigenvectors. The same is true of 
the nearest- neighbor pieces of A ia jp, by far its largest 
matrix elements, which are also proportional to J . The 
net result is that the lower pair of flat bands remain flat 
while the upper pair acquire the same dispersion as J — 
this is illustrated in the middle panel of Fig. 

Thus far we have shown that "much" of H is charac- 
terized by a pair of low lying flat bands. What remain, 
the matrix elements of A ia jp beyond the nearest-neighbor 
distance, are small. Thus their inclusion does not mod- 
ify the gross features significantly, preserving a spectrum 
similar to that of the nearest-neighbor model. Indeed, 
adding Ay 8 beyond the nearest-neighbor distance weakly 
splits the remaining pair of flat bands, which acquire a 
small dispersion; the other pair of bands is barely modi- 
fied on the scale of their dispersion; this is shown in the 
figure's bottom panel. 

The full A was included by adding its numerical 
Fourier transform to that of J and V obtained analyti- 
cally, and diagonalising the resulting 4x4 matrix. With 
Ai Q j/3 small and decaying fast, its Fourier transform is 
quickly and absolutely convergent; we have checked that 
the results are essentially independent of the truncation 
distance for r c > 12r nn . We would like to stress that 
this truncation is not equivalent to the truncation of the 
long-range dipole interaction since in our case the main 
long distance part of the dipole interaction is already 
contained in V. 

In Fig. |31 we show the band of maximal eigenvalues of 
the dipole matrix. It is weakly dispersive because A^ 8 
is small. Our results are in quantitative agreement with 
the Ewald summation results [n] • 

Ground states and entropy: The nearest neighbor 
interaction by itself gives rise to a macroscopic entropy at 
T = 0. We now show that any combination of the nearest 
neighbour, J and the model dipole, V, interactions leads 
to precisely the same ground states, and hence entropy. 

First, note that the ground states of J have zero en- 
ergy; an arbitrary ground state, |0) can thus be written 
as a linear combination of modes in the flat zero-energy 
bands only, as an admixture of other bands would lead 
to a non-zero energy: 

2 

l°> = EEMq)Mq)> ■ (5) 

q fj,= l 

The hard-spin condition imposes constraints on the am- 
plitudes a M (q), which we do not need to resolve explicitly. 

It suffices to note that, having written the ground 
states of J in terms of the flat bands, it follows that they 
remain ground states upon admixing the model dipole in- 
teraction, V , as this only affects the relative position of 
the excluded positive energy bands /i = 3,4. Indeed, 
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FIG. 2: (color online). Mode spectrum of the dipole inter- 
action (the eigenvalues, Vi, of the Dr^ n T>fP matrix). Top: 
Model dipole interaction, V only. There are two degener- 
ate pairs of flat bands. Middle: After inclusion of nearest- 
neighbor correction. Two bands remain degenerate and flat. 
The other two become dispersive. Adding a nearest-neighbor 
superexchange can enhance/suppress this change. Bottom: 
Correlation function plus A^ 3 . The remaining flat bands be- 
come weakly dispersive. The middle and bottom plots are 
almost identical on this scale. (The spikes in the figure have 
to with the special character of the single point at q = for 
the conditionally convergent dipole interaction.) 

this implies that the set of ground states of J and V are 
identical. 

The addition of A will lift this degeneracy but, as we 
saw, only weakly. Consequently the low temperature en- 
tropy of spin ice will be quite close to that of ice before 
it goes away at the lowest temperatures as the system 
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FIG. 3: (color online). Maximum eigenvalue, v max , of the 
dipole matrix in the [hhl] plane, q is in units of 2-k, r c = 16r nn . 

orders 0. 

Correlations: This explains why the similarity of 
the spectrum is sufficient to yield the correct low- 
temperature physics: for the model dipole interaction, 
the ground state correlations are exactly those averaged 
over the ice rule manifold. However, this equivalence be- 
tween the model dipole and the nearest neighbor antifer- 
romagnetic Ising model breaks down at nonzero tempera- 
tures. For the nearest-neighbor Ising model, the presence 
of thermally activated ice rule violating defects leads to 
an exponential decay of correlations on a length scale 
diverging as £ ~ exp(2J/3T) at low T. 

By contrast, for the dipole problem the long range of 
the interaction implies long ranged correlations at any 
temperature. This is already evidenced by the first term 
in the high-temperature expansion of 

(S ia ■ 8 jjfl ) cx -H ia , jf3 /T ~ O(r-r^) . ( 6 ) 

In fact, in a saddle-point treatment for V , one can show 
that this holds for any temperature, and to all orders in 
corrections to the saddle point. 

Projective equivalence: The mathematics under- 
lying our analysis of dipolar spin ice can be general- 
ized. One can construct other exchange matrices J' 
which share their low-lying flat bands (and its eigenvec- 
tors) with interactions V of longer range; we should note 
though, that generically, neither the J 1 nor the V are 
of bounded range. In this fashion, we find pairs of inter- 
actions which are equivalent under projection to the flat 
bands; this we term "projectively equivalent" . 

The miracle of spin ice is hence twofold: first, that 
the physical dipole interaction restricted to the site de- 
pendent easy axes on the pyrochlore lattice provides one 
member of such a pair; and second, thanks to its emer- 
gent gauge structure, the other member of the pair is the 
classic ice problem dating back to Bernal, Fowler and 
Pauling. In short: dipolar spins are ice because ice is 
dipolar. 
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